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Abstract. Robustness of the solutions to the inflaton potential inverse problem based 
on the slow-roll approximation is addressed. With that aim it is introduced a measure 
of the difference of the outputs obtained using first and second order respectively in 
the horizon-flow expansion. The evolution of this measure is determined by a second 
order linear non-autonomous non-homogeneous differential equation. Boundedness of 
the general solutions to this equation is analyzed. It is shown that they diverge for 
most of the physically meaningful cases. Examples for typical inflationary models are 
presented. It is argued that this lack of robustness is due to the limitations of the 
slow-roll expansion for probing the scale-dependence of the inflationary spectra. 
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1. Introduction 

The relatively recent release of the analysis of the first year WMAP data [Tj and of a 
significant amount of very precise data of the power spectrum of matter distribution 
at large scales j2j caused a turmoil in the cosmological community. The feasibility 
of obtaining reliable information about our universe when it was younger than 10~^° 
seconds accounts for a large part of that excitement. Data confirmed that the primordial 
power spectrum of curvature fluctuations is consistent with a scale-invariant, Gaussian 
and adiabatic spectrum. The simplest and most elegant mechanism known to causally 
produce primordial spectra with those properties is cosmological inflation. 

In the simplest scenarios with an exit from inflation the responsible for the 
accelerated expansion of the early Universe is the inflaton ip, a single scalar field with 
equations of motion, 

ip + 3Hip + V'{^) = 0, (1) 

(2) 



Here V{if) is the potential energy density of the inflaton field, dot stands for a derivative 
with respect to cosmic time t, and prime denotes a derivative with respect to ip. We use 
Planck units, SttG = c = h = 1. 

The Hubble horizon c/h = ^/H, is roughly the size of the region where causal 
processes can take place during one Hubble time 1/H. During inflation the comoving c^h 
decreased allowing comoving scales to cross out the causal horizon. Physical information 
of the inflaton and of the space-time quantum fluctuations contained on those scales 
was effectively 'frozen'. This way, when the scales crossed back later during the era 
of standard (non-accelerated) expansion, the whole observable universe was seed with 
statistically uniform curvature perturbations for large-scale structure formation. These 
seeds served also as initial conditions for the evolution of CMBR anisotropics. 

The spectra of quantum fluctuations of density and metrics during cosmological 
inflation can be predicted [3]. It is required to know the behavior of the Hubble horizon 
as a function of the logarithm of the scale factor, N = Ina(t). In the inflaton scenario 
it means solving Eqs.(0) and (j2I). For most potentials that is a difficult task that can 
be simplified using the horizon-fiow functions which are defined recursively as jl]: 

Wi = ^^ Vm>0, eo = ^, (3) 

where A^i denotes an arbitrary 'initial' moment. The necessary condition for infiation 
to take place is ei < 1, and if the weak energy condition and null energy condition hold 
true, ei > 0. In general, for the primordial spectra to be nearly scale invariant, it is 
sufficient to assume \em\ < 1 allowing, this way, to expand all the interesting infiationary 
quantities in terms of the horizon- flow functions^. The slow-roll approximation is a 
particular case where all the are assumed to have the same order of magnitude. 

I Note that according to recent results that condition may be unnecessary |S] . 
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Eqs.(PQ) and relate the horizon-flow functions to the inflaton potential V{(p) and 
its derivatives with respect to the inflaton fleld. For instance, the flrst two flow functions 
are given approximately by 



Thus, approximated expressions for the inflationary spectra can be derived. To flrst 
order for the horizon-flow functions§ they were found by Stewart and Lyth [7j using 
the slow-roll approximation to truncate the argument of Bessel functions evaluated at 
the time of horizon-crossing. Their accuracy was tested against numerical solutions in 
Refs. jSl 13 El- In the same approximation but using Green's functions, second order 
expressions were derived in Refs. [101 Ej and tested by Schwarz and Terrero jSj . The 
approach based on Bessel functions was then modifled jH IHI to derive expressions at any 
order for two wide classes of inflationary models that do not necessarily satisfy the slow- 
roll condition. Up to third order, the corresponding expressions were tested in Ref. [Hj. 
Several other methods were introduced to allow for relaxing the standard slow-roll 
condition fTH, IT^ ITH] . In general, their predictions seem to match the accuracy of those 
methods tested against numerical results. All these tests yielded that the predictions of 
inflationary spectra using terms up to second order in the horizon-flow expansion are 
accurate enough to match the precision of current and near future observations||. 

This way, given an inflaton potential as input in Eqs.(^ and (j21), the horizon- flow 
functions (jH)) can be evaluated at any moment, the behavior of the Hubble radius can be 
described by means of (0), the spectra of the inflaton and metrics fluctuations accurately 
predicted, used as initial conditions for the calculation of the evolution of scalar and 
tensor modes of anisotropies of the CMBR, and the resulting spectra compared with 
the observational data to test the reliability of the given inflationary model. 

'Observing' the inflaton potential means solving the related inverse problem, i.e., 
that of looking for the functional form of the inflaton potential corresponding to given 
CMBR spectra. Herein, we will refer to this task as the Inverse Problem for the Inflaton 
Potential (IP^). Starting with the seminal paper by Hodges and Blumenthal ^1], during 
the last 15 years there were introduced several methods for solving the IP^. Depending 
on the way the output is given, they can be classifled in parametric [T ^ ITK | \W \ ITT f IT^. 
perturbative [1303120], full-numerical [21], and stochastic [22] methods. 

Assuming that the initial conditions for the calculation of CMBR spectra have been 
already fltted to the observational data, what all the methods for solving the IP^ do (with 
the remarkable exception of the full-numerical method) is, essentially, to start by looking 

§ Since the order is not the same for aU the observables, we follow here the convention in Ref. [U] of 
denoting the order by that of the first term in the Taylor parametrization of the primordial spectra. 
II There are a few issues on this regard that will be commented in the last section of this paper. 





In turn, the behavior of d}i{N) near a suitable pivot point N^^ is described by 
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for the combinations of that gives the best fit to those primordial spectra. These 
combinations are constrained by the form and order of the approximated expressions for 
the spectra of density and metrics fluctuations. Note that, following definitions (jH)), the 
number m of horizon-flow functions is directly related to the order of the corresponding 
expressions. In its turn, the order of the predicted spectra is constrained by the finite 
precision of data. Now, according with expansion (0), specifying the set {em{N^)} is 
equivalent to specifying d}i{N). So, the unavoidable truncation of {em{N^)} corresponds 
to an incomplete knowledge of the evolution of the Hubble horizon. The next step in the 
IP^ solution is to use Eqs.(IT)) and (0), or the approximations (jH), to find the functional 
form of the inflaton potential corresponding to that approximated behavior of dniN). 

The reliability of solving the IP^ for a potential as constrained by measurement 
errors was directly tested in Refs. [201 EHj. There is also a large number of studies on 
the impact of the observational uncertainty on the form of the output potential. They 
give us complementary information about possible troubles for observing the potential 
(see Refs. j21l El El 123 ESI EZI and references therein). Though the results are in general 
encouraging, it must be noted that each of them uses in their analysis a fixed order for 
the slow-roll expansion of the primordial spectra. This is a formal expansion, hence there 
is not rigorous proof that it converges so, there is not certainty about the solution of the 
inverse problem based on this expansion to be robust. By robust we mean that, while 
increasing the order of the underlying expressions, the IP^ solution must converge to a 
unique functional form. This robustness is required in order to make definite statements 
about the high energy physics linked to the observed inflaton potential. 

In this sense, the preliminary analysis presented in Ref. [2H1 gave us a first 
cautionary warning about the possible existence of troubles for the robust solution of the 
IP^ when the slow-roll approximation is used. The aim of this paper is to step forward 
in the analysis of the robustness of the reconstruction of the inflaton potential with 
regards to the convergence of the horizon-flow expansion for the spectra of inflationary 
perturbations under the slow-roll approximation. With that in mind, in the next section 
we will derive the basic equation for our analysis, a second order linear non-autonomous 
non- homogeneous differential equation obtained from the comparison between equations 
for the inflationary perturbations to first and second order in the horizon-flow expansion. 
In SecEjwe will present mathematical evidence pointing to the non-existence of bounded 
solutions for this equation if the time variation of the coefficients is required to be 
consistent with the inflationary cosmology. We reinforce that conclusion in SecEl based 
on the analysis of inflationary models belonging to the prototypical classes introduced 
in Ref. [Oj. Finally, in Sec|S]we discuss our results. 

2. The basic equation 

For finding the functional form of the infiaton potential all the methods mentioned in the 
Introduction use the same kind of expressions for the spectra of scalar perturbations. 
Nevertheless, early in the research on the IP^ it was realized that the role of tensor 
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perturbations deserves also attention if a unique solution to that problem is desired j2Sl- 
The information on these modes in the available data is so far rather poor ^1 12] . That is 
why all the studies but the so-called Stewart-Lyth inverse problem [THj avoided to deal 
with the tensor equations. However, the quality of the data on the tensor modes could be 
radically improved if a number of interesting experiments manages to see the primordial 
light (see for instances Refs. [221 )• Even more, studies based on the parametric method 
have shown that even a poor information on the spectrum of the primordial gravitational 
waves allows for breaking the degeneracy in the solution of the IP^ and for obtaining 
interesting results [IHl HH IS2] • Particularly, in Ref . ^| it was shown how using the 
ratio of the amplitudes of the tensor and scalar perturbations, 

r = a—^ , (tt = constant) , 

as input for the IP^ can be a very fruitful way of taking into account the available 
observational information on both types of spectra to observe the convexity of the 
potential during inflation. To do that it is necessary to solve the differential equation 
(recall the definitions Q for the e^). 

In — = In ei + Ce2 + CidCs + C2el + C3e2e3 , (6) 

where rg is a constant and 

C ^ -0.7296 , (7) 

Ci = l^-y + 5 + ^ -0.6644 , (8) 
C2=(^-^ + l^ ^-0.2337, (9) 

This equation can be derived from the second order expressions for the amplitudes of 
the scalar and tensor perturbations obtained using the slow-roll approximation |10[ lUj. 
(It can also be derived following the lines described in Ref. 

If the reconstruction of the infiaton potential were robust, then we would expect 
that, given r/ro as input into the IP^, the output using up to first order in Eq. © (see 
Ref. jn2| for an example of such an output) will differ only in small features from the 
output of the problem with the same input but taking into account all terms in To 
assess this small departure, following Ref. [2^, we introduce a quantity S{N) such that 

%i = ei(l + 5), (11) 
where the 's' denotes the horizon-flow functions in the set to be found if second order 
terms in Eq. © were used. The standard notation remains for the set found with up 
to first order terms. Next, using definitions (jH)) it can be obtained 

(12) 

%2%3 = e2e3 + ^, (13) 
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Now, by definition ei, €2 and €3 satisfy identically the first order version of Eq. ^ 
and we will demand S{N) to remain very small as increases. Substituting (fTT|) . (fT^ 
and (|T!?|l into (jHl), and keeping only linear terms in 6{N) and its derivatives, the following 
second order linear non-autonomous non- homogeneous differential equation |28j, 

^.«(.,^-.^.(.)^/(.,, 

is obtained, where 

„(^)^(£±£i^,±^, (15, 

cu^ = --^ ^ 6.893, (16) 

/(AT) = (Cieie2 + C2el + ^36263) ^^^^ 

C3 

For the inverse problem of the infiaton potential to be robust, the mildest 
requirement we could ask for the solution of this equation is to be bounded for all 
possible initial conditions S{N[). As we will show in the next two sections, that seems 
to be unlikely for all a{N) and f{N) physically meaningful. 

3. Mathematical evidence 

Equations like ()14|1 are usually found in mechanical problems like the study of vibrations. 
As in Mechanics, a{N) is going to be called here the 'damping' of the system, —u"^ the 
'stiffness', while f{N) will be referred as the 'forcing'. 

First of all, we recall that the formal solution of Eq. ()14|) is 

N 



6iN)=g''5iN;)+ igT'fis)ds. 
where : is the linear (A^i, A^)-advance mapping for the homogeneous system 

As we already mentioned, for most of the cosmologically interesting cases, \em\ < 1 for 
m > 0. According with definitions Q this means that they vary slowly with A^. So, 
without lost of generality, we can assume the forcing f{N) to be bounded. In that 
case, boundedness of ffTHj) requires boundedness of solutions to Eq. (fT^. The necessary 
conditions for the existence of bounded solutions of this equation with undetermined 
coefficients is still an open problem. Nevertheless, there are inferences that can be drawn 
for the case of a stiffness strictly negative for all A^. Such a stiffness is usually considered 
as a signal of instability of the trivial solution of the above equation. The reason is that 
for the corresponding to Eq. (fT^ planar system, 

1 I 

' A with A^ = \6{N) 



dA 



dN 



a;2 - a{N) 



T _ i./^n d6{N) 



dN 
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the eigenvalues A^(A^) oc y—a{N) ± a{Ny + Au'^j point to a saddle-like instability 
for any given value of A^. If these eigenvalues are slowly changing, then the system 
is assumed to have no bounded (finite) solutions. A rigorous foundation of this folk 
theorem is still missing together with a precise definition of 'slowly changing'. However, 
Lyapunov's direct method reveals that, for the existence of bounded solutions to 
equation |T9|l for all S{Ni), the following conditions are close to necessary jHlj: the 
damping a{N) should be positive (at least asymptotically) and should increase in an 
exponential rate. As it was already mentioned, this last condition is incompatible with 
all physically meaningful inflationary scenarios. 



4. Physical evidence 

To consolidate the analysis in the previous section, let us see how the solutions of 
equation (fT^ behave for physical situations. To cover as much as possible the space 
of inflationary scenarios we will use one model from each of the classes introduced in 
Ref. p. We choose that classiflcation because it involves only exact expressions. Then, 
the classifying criteria are not dependent on the order of the horizon-flow expansion, 
contrary to what happens in the more frequently used small- fleld/large- field/hybrid 
models classiflcation [25- This was realized flrst by Schwarz and Terrero and recently 
conflrmed by Kinney and Riotto [2Z|. Classiflcation in Ref. [5] is unambiguous because 
is strictly based on physical criteria, namely the behavior of the kinetic and total energy 
densities of the inflaton fleld. It consists of three classes, I) hidden-exit inflation (e2 < 0), 
II) toward-exit inflation with general initial conditions (0 < e2 < 2ei), and III) toward- 
exit inflation with special initial conditions (0 < 2ei < €2). 
Monomial potentials like 

V(^) = A— (20) 
n 

with chaotic initial conditions give rise to models of II kind [35j. Using relations (ji} it 
is found that for such potentials 

ei(v^)-y^. 

To convert from dependence on ip into dependence on the efolds number, we use 

rJrh , 

±v^, (21) 



dN 

what can be derived using deflnition (jS)) for ei and the Friedmann equation ^ for the 

inflaton cosmology. After doing the conversion it is obtained, 

n 4 
ei ^ - {Nf -N) , 62 = 63 = -ei . 
4 n 

So, the class criterion is met for n > 2. Nf is the number of efolds by the end of inflation. 
Substituting the above results into (|TT)|l and (|T7jl the following expressions are obtained, 

nCi + 8C2 C , , 
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fiN) 



nCi + 4C2 + 4C3 



(23) 



ACsiNf-Ny ■ 

Next, we derive the analogous expressions for class III (toward-exit inflation with 
special initial conditions). Typical for this class are models where the inflationary phase 
starts near a false vacuum. An example arising in Superstring Theory is 



V = Vo- 



Here, according to ( 



in the slow- roll phase we have. 



(Vo 



2 ' 



(-2 



ei . 



This kind of models exits inflation gracefully when reaching ei{Nf) = 1 so, for simplicity, 
we set the condition (fiNf) = ^/2Vo/m? (implying Vq < 1) and use (|7H) to obtain. 



ei ~ exp 



2m^ 
1^ 



{Nf - N) 



(-2 



2m^ 



63 = 0. 



Then, the corresponding damping ()16|) and forcing ()17|) are given by. 



a{N) 



fiN) 



9i 



■ exp 



2m^ 



exp 



(Nf - N) 



AC2 2m2 C 



2m^ 
1^ 



(Nf - N) 



AC2m^ 



(24) 



(25) 



Finally, typical examples from class I (hidden-exit inflation) are hybrid models 
where cosmic acceleration is driven by the inflaton and, in order to end inflation, there 
is yet another field which finally drives the vacuum energy Vq to zero. Our example is 
the hybrid model with effective potential [HZ| 



V = Vo + 



2 ' 



€2 



Vo 



In this case from we obtain 

1 m^(j? 

2 (Vo + 

Similarly to the previous example, but now in the opposite regime, to simplify the 
calculations we assume the initial conditions for the inflationary era^ ^i{N{) = 1, 
^p{Ni) = v^Vo/m^, what after using (pT|l leads to 



ei ~ exp 



2m^ 
1^ 



(iv - ivo 



£2 



2m^ 
1^ 



63 = 0. 



The corresponding damping and forcing are then given by. 



a{N) 



fiN) 



■ exp 



C3 
2Cim^ 



2m^ 



Vo 



{N - Ni) 



4^2 2m2 C 



exp 



2m^ 
"1^ 



{N - N,) 



(26) 



(27) 



% There are uncountable others choices that will led to similar results and we should expect our analysis 
to be independent on the initial conditions for inflation. 
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Figure 1. Plot of the solutions of equation l|14() for the chaotic model V{(p) = Xip^/A. 
Each curve corresponds to a different solution starting at A^i = (blue) and running 
towards Nf = 100 (red). 



After substituting expressions and (|^. and putting the numbers (ITUj) . one 
can look for numerical solutions to Eq. In Figgis shown the phase-space for a 

potential given by Eq. (j^n|) with n = A. Note that the value of n does not modify the 
qualitative behavior of these solutions, while the value of A does not play any role at 
all. Hence, we can conclude that for this kind of models the solution obtained to second 
order in the horizon-flow functions will typically depart from the lowest orders solutions. 

The phase diagrams for the damping and forcing given correspondingly by (j24j) and 
(123), or by and ^T7\ . are topologically equivalent to that in figure [TJ Therefore, 
it is confirmed the result in the previous section about solutions to Eq. (fT^ diverging 
around a saddle-point near the origin, as the generic dynamical pattern for those cases 
which are interesting from the point of view of inflationary cosmology. 



5. Discussion 



The rigorous conclusion to be drawn from the results presented in this paper is that 
solutions of the inverse problem for the infiaton potential obtained using expressions up 
to second order in the slow-roll approximation will be generically very different from 
those solutions obtained using lower order expressions. 

It may happen that including higher orders helps to get rid of this trouble. However, 
the order of the expressions used to predict the inflationary spectra is determined by the 
accuracy of the data from the observations of the CMBR anisotropies and large-scale 
matter distribution. As it has been pointed out, second order expressions will indeed 
match the precision of current and near future measurements jHlEl- In order words. 
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in the expressions for the amphtudes of the inflationary spectra, horizon-flow functions 
with orders higher than 2 would not contain any useful information. Therefore, including 
such orders would be irrelevant for the output of the inverse problem. 

Note now that equation (fT^ has trivial solution S{N) = only for the cases of 
de Sitter (where = for m > ) and power-law inflation jHEI- In this last model 
a (X (with a constant p ^ 1), implying oc exp(iV/p), ei = 1/p, = for 
m > 1 and V{lp) = Aexp[B(p). Here B is a. real constant which is zero for a universe 
dominated by a cosmological constant. Therefore, S{N) actually measures the depart 
from an exponential of the functional form of the output potential. Power-law inflation 
is the unique model which predicts exact power laws for the amplitudes of both the 
tensor and scalar spectra i.e., A{k) = A^.k"', where k is the comoving wavenumber 
and n is a constant called the spectral index. Moreover, during power-law inflation 
the scalar index ns is equal to the tensor ut, what implies a constant ratio of tensor 
to scalar amplitudes, r. Conversely, in Ref. [201 it was proved that using power- law 
spectra as input in the inverse problem renders the exponential potential as the unique 
output. From this point of view, S{N) measures the depart of the behavior of the 
inflationary spectra from a power-law. We can conclude that the essence of the troubles 
for observing the inflaton potential is the weak signal in the CMBR of scale-dependence 
of the primordial spectra and the difficulties to extract that information from data using 
expressions for the inflationary perturbations based on the slow-roll approximation. 

As a matter of fact, we were able to clarify here a point which evidences can be 
traced back in the existing research literature. In 1995, Lidsey et al. [201 showed that 
the results for the spectra obtained to first order in the slow-roll approximation by 
Stewart and Lyth [2j could be understood as an expansion about the exact solution for 
power-law inflation. Two years later, Wang et al. [IHl analyzed the implications that 
finite orders in such an expansion have for precisely predicting the inflationary spectra. 
Limitations pointed out were linked to the Bessel and horizon-crossing approximations. 
Martin et al. jH] quantifled these limitations by comparing a power-law parametrization 
for the amplitudes with a more general Taylor series parametrization, and by analyzing 
the impact of the variation of the pivot scale on the accuracy of predictions. The larger 
the scale-dependence of the spectra, the bigger is this impact (see also Ref. [27]). 

The combination of the Bessel and slow-roll approximations is no longer needed to 
calculate the primordial spectra, and the power-law parametrization has been extended 
to included the running of the scalar index. However, the slow-roll condition remains on 
the basis of the expressions that usually link the Taylor parametrized spectra and the 
inflationary predictions. There is evidence that the above mentioned limitations still 
persist. Comparing approximated and numerical predictions for three given models, in 
Ref. [H] it was noted that increasing the order of the horizon-flow functions beyond the 
quadratic in every included term of the Taylor parametrization of primordial spectra 
does not signiflcantly improve the accuracy of the approximations. In fact, adding terms 
in that parametrization actually seems to be most important for increasing the precision 
of the predictions. However, the number of terms in the parametrization depends 
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on the capability of the observations to measure a significant departure from scale- 
invariance in the primordial spectra. As a final evidence, we recall that using first order 
expressions for r and ns, Dodelson et al. were able to discriminate among inflationary 
models with decreasing and increasing kinetic energy density by distributing them in 
the (r, ns) plane [24J. As it was noted in Ref. jH], including higher order corrections 
makes impossible to distinguish between these models on the basis of a r — ns plot, 
unless the running of the spectral index is known or constrained to be very small. The 
same conclusion was recently made by Kinney and Riotto |27j . 

Limitations from theoretical and observational sides to account for the scale- 
dependence of the primordial spectra have been misunderstood as generic predictions 
of inflation. That this scale-dependence is in the edge of observational capabilities does 
not means that all the successful infiaton potentials must resemblance an exponential 
during the inflationary era, but just that we are unable so far to observe the difference. 
The findings by Hoffman and Turner jl^I about power-law inflation as an a-attractor 
(attractor in the future) of the flow in the space of inflationary models were constrained 
by the order of the expressions in their analysis. As it was shown in Ref. , departure 
from power-law behavior is nullified whenever the leading order is used. To first order, 
for instance, it remains true that the inflationary dynamics yielding a constant tensor 
index nx has power-law inflation as an a-attractor jSO] . Nevertheless, it is also true that 
the inflationary dynamics characterized by a constant ratio r, has power- law inflation 
as an a;-attractor (attractor in the past or repellor) Furthermore, in Ref. it 

was shown how a minimal scale-dependence for r (and information on such a scale- 
dependence can be drawn from the difference between both spectral indices) allows 
for models with power-law inflation as a transient regime. Based on these facts, we 
conjecture here that, if using higher orders in the slow-roll expansion for the analysis of 
the inflationary flow, power-law inflation still arises as a fixed point, then it will do it 
like a saddle-point, consistently with the results presented in this paper. 

Summarizing, either the slow-roll expansion does not converge or it does so 
slowly that a large number of terms in this expansion is required to probe the scale 
dependence of the primordial spectra. Nevertheless, information on this dependence 
seems to be a necessary ingredient for a robust solution to the infiaton potential inverse 
problem. This conclusion should encourage further development of the full-numerical 
method for the resolution of this inverse problem [21], as well as of other analytical 
methods, for instance, those based on the WKB ^T] approximation, on the uniform 
approximation ^2), and on the inverse-scattering theory 
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